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Abstract 

We present strong evidence that a coupled-map-lattice model for spatio- 
temporal intermittency belongs to the universality class of directed percola¬ 
tion when the updating rules are asynchronous, i.e. when only one randomly 
chosen site is evolved at each time step. In contrast, when the system is sub¬ 
jected to parallel updating, available numerical evidence suggest that it does 
not belong to this universality class and that it is not even universal. We 
argue that in the absence of periodic external forcing, the asynchronous rule 
is the more physical. 
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The onset of spatio-temporal intermittency is a common phenomenon of many extended 
systems ranging from models based on conpled-map-lattices to varions experiments 
in convection [^|[ and in the “printers instability” p,|^. A particnlar elegant conpled- 
map-lattice (CML) showing spatio-temporal intermittency was introdnced some years ago 
by Chate and Manneville This CML employs individnal maps that can be in two 

very different states: either in a chaotic (or “tnrbnlent”) state or in a “laminar” state. For a 
single map the laminar state is “absorbing”: once the motion is in the laminar state it cannot 
escape. For the conpled system, one observes interesting dynamical patterns of tnrbnlent 
patches penetrating into a laminar state, and becanse of the strong flnctnations, this has 
been called |]l| spatio-temporal intermittency (STI). Once the system is in a pnre laminar 
state, it cannot escape: this is an absorbing state of the fnll spatially conpled system. These 
properties of the STI led Pomean [§ to conjectnre that the critical properties at the onset of 
STI should be governed by the exponents of directed percolation. The turbulent spots of the 
dynamics in a space-time plot percolate through the system in a manner very similar to the 
connected bonds in directed percolation (DP), which also has an absorbing state. Subsequent 
extensive numerical studies and scaling arguments did not show agreement with this 

conjecture. On the contrary, it was found that the generic critical properties were not in 
the universality class of directed percolation. In fact, since the critical properties vary with 
the parameters, they are not even universal. Nevertheless it was argued in that the 
apparent non-universality is due to travelling solitary excitations with long life times, and 
that one should in principle recover the DP-behaviour only for extremely long time scales. 

The standard time evolution of a coupled-map-lattice is by synchronous (or parallel) 
updating, in which all individual maps of the lattice are iterated forward simultaneously 
in a completely deterministic way The connection between such coupled map lattices 
and physical systems, described by partial differential equations is, however, not straight 
forward. In particular, the synchronous update can only be motivated when the system 
dynamics is driven by a global periodic external force (or “clock”). Recently it was observed 
that the critical properties of a standard CML, using coupled logistic maps, actually depend 
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on whether synchronous or asynchronous updating is applied |T^. Here, asynchronous means 
that in each step a random site on the CML is chosen and is iterated forward. 


This result forces us to reconsider the evidence for non-universal behaviour for STI, which 
is all based on the synchronous updating rule. In this letter, we thus apply the asynchronous 
updating to the coupled-map-lattice of spatio-temporal intermittency discussed above. We 
hnd that all critical exponents, independent of the choice of parameters, fall into the univer¬ 
sality class of directed percolation. All critical characteristics of DP, such as hyper-scaling, 
seem to be fulfilled, leading to independent controls of the values of the critical exponents. 

The original dynamics of the coupled-map-lattice of Chate and Manneville |]I| for a lattice 
with one space- and one time-direction is written in terms of a field Ui{t) at site i and time 
t as 


Ui{t + 1 ) = + ^AfUi{t ), 


( 1 ) 


where AfUi(t) = f{ui-i(t)) —2f{ui(t)) + f{ui+i(t)). The parameter e measures the strength 
of the coupling from site i to its two neighbours. The dynamics (2) is parallel or synchronous: 
all sites are updated at the same time. 

The local map / is of the form: 


rx, if a; G [0,1/2]; 

/(x) = <1 r(l-x), ifxe[l/2,1]; (2) 

^x, if a: g] 1, r/2]. 

The chaotic motion of / for a; < 1 is governed by a standard tent map of slope r. However, 
when r exceeds the value two, the trajectory may escape to a “laminar” state with x > 1, 
and this state is marginally stable, because the slope in the line of fixed points is one. 
As mentioned above, the laminar state is absorbing, i.e. the trajectory cannot be pulled 
back into the chaotic state. This is no longer the case, when the maps are coupled, since 
the interactions with its neighbours may pull a laminar site back into chaotic motion, thus 
causing the interesting interplay between laminar and turbulent regions. 

In this work we consider the same CML under asynchronous updating: at time t a 
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random site v is selected and is altered according to (|I]) while all the other variables keep 
their values: 

Ui^{t + 1/L) = f{ui^{t)) + , (3) 

Ui{t + 1 /L) = Ui{t) for i ^ ir , 


where L is the size of the system. Note that with this choice of time step on average each 
site is updated once in one unit of time. 

The introduction of chaotic iterated maps in extended dynamical systems - and especially 
simple one-dimensional, non-invertible maps - can only be motivated in a very heuristic way 
|jl|]. The closest physical parallel is probably a collection of weakly coupled subunits, each 
able to perform chaotic dynamics. In order to derive a discrete, local map, one applies 
a Poincare section. Note that the Poincare map owes its simplicity to the fact that it 
is obtained by restricting the dynamics to a surface through the local phase space and 
is therefore not in general equivalent to moving the system forward through a fixed time 
interval. In the absence of an external, periodic forcing, the time interval between consecutive 
crossings of the section will generically be variable and thus vary from unit to unit in space. 
To find the state of the entire system at fixed time intervals, this variation has to be taken 
into account. This can, in a rough way, be accomplished by using the asynchronous update 
in which the units experience slightly different update times. Of course we are thereby 
replacing the complicated, deterministic update rule by a random one and this is obviously 
a crude approximation. But we believe that this is often closer to physical reality than the 
synchronous update. As an example, synchronous updating can lead to complicated, spatial 
structures, in cellular automata [|1^, which disappear under asynchronous updating |1l4| , 
and should therefore be considered unphysical in the absence an external clock. 

Fig. g shows a pattern generated by the asynchronous update introduced above. As 
only one site is updated at a time, a new horizontal line in the time axis is added only after 
L = 128 time steps (i.e. t ^ t + 1). We observe that the turbulent sites (black in Fig. 
percolate through the system, sometimes ending in a dangling bond. 
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To estimate the critical exponents for the randomly updated CML (^, we follow the 
hnite size scaling methods of Houlrik et ah P, p!5| , p!6| . First of all we have to locate the 
critical line in the parameter plane (e,r). This is done by measuring the absorption time 
r(r, e,L), i.e. the time it takes the system starting from a random initial state to reach 
the absorbing state - averaged over an ensemble of initial conditions. At the critical point 
e = ec(r), this time diverges like 

r(e„L)~L^ (4) 

where the usual dynamical exponent z = has been introduced. Fig. || shows the phase 

diagram with the critical line and contrasts it with the critical line for the synchronously 
updated system ([^) taken from |^. We consider three different values of r and the corre¬ 
sponding values of Cc and 2 ; are found in Table |. 

The order parameter m(e, L, t) is dehned as the fraction of turbulent sites in the lattice, 
again averaged over many different initial states [Q. The order parameter scales in the usual 
way when approaching the critical line from above 

m ~ (e — ecY for ^ ^ (5) 

Using this relation we have estimated the /3-exponent directly and the results are found in 
Table |. Applying hnite-size scaling arguments we hnd the following scaling form at the 
critical point 

The function g{t/L^) has the following properties: At times much smaller than the ab¬ 
sorption time, one expects a power-law decay in time due to critical correlations and the 
L-dependent pre-factor in (j^) must drop out. For t L^, we therefore have 

m(ec,f,T) ~ (7) 

For times much larger that the absorption time, we may assume uncorrelated decay of the 
order parameter, and therefore the function g, will decay exponentially. Fig. shows a plot 
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of m{ec,t,L) versus t at r = 2.2. The curves for 7 different system sizes in the interval 
L = 2^,...., 2^° fall very accurately on the same line in the double logarithmic plot allowing 
determination of as listed in Table |. As we see in this table all the critical exponents 
are, within the error bars, consistent with the values for directed percolation, listed at the 
bottom line in Table |. 

In order to obtain independent checks on the values of the critical exponents we have 
performed a rescaling analysis using (|^). The rescaled curves are shown in Fig. For 

systems sizes in the interval L = 2"^,...., 2^° the rescaled curves collapse very accurately to a 
single curve. The corresponding values oi j3/v± and 2 ; are shown in Table 0- again consistent 
with DP. 

To extract further critical exponents and to check hyper-scaling we have looked at the 
spatial correlations. Even though any dynamics of this CML will end up in the absorbing 
(laminar) state, it is possible to hnd non-trivial spatial correlations in a long lasting quasi 
stationary situation []^. The correlations can be obtained from the pair correlation function 


- («W)^ 

^ i=i 


( 8 ) 


where the brackets denote the average over different initial conditions. If the coupling 
between the sites is weak, i.e. if e is small, one might expect the spatial correlations to fall 
off exponentially with distance. At criticality, on the other hand, one expects an algebraic 


decay of correlations |15 


Cj{t) =j^ , 


(9) 


where rj is the associated critical exponent, and correlations are induced over a length scale 

m fl/z 

as the CML relaxes towards a steady state. The (equal time) correlation function 
for various times (at r = 2.2) is shown in Fig. indicating an algebraic decay in space after 
long time. The corresponding value of r] from this direct measurements is shown in Table | 
together with values obtained at other r-values. Again the agreement with DP is conhrmed. 
For the spatial correlations, a rescaling analysis can also be performed by plotting j /A 
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versus Fig. shows the corresponding rescaled plot and this technique allows 

another independent determination of z and r], the values of which are shown in Table 
The fourth column of Table contains the values of rj obtained from the hyper-scaling 
relation 


2(5/u = d-2 + r] , (10) 

giving a third way of estimating the exponent rj. All three ways of hnding rj give results 
that, within error, are in agreement with the directed percolation value 1.51. 

We thus conclude that our numerics gives very strong evidence for the fact that spatio- 
temporal intermittency in the form of asynchronous coupled maps falls in the universality 
class of directed percolation. The reason why synchronously updated maps do not behave 
in a universal way must thus be sought in the complicated correlations built up by the 
strong constraint of exactly simultaneous, completely deterministic updating, which, in most 
applications will not be fulhlled. Thus, experiments showing spatio-temporal intermittency 
should be describable by directed percolation near the transition to turbulence - at least 
insofar as they do not involve global periodic forcing. One might speculate on the possibility 
of observing the non-universal “synchronous” behaviour in periodically forced systems, as 
say the printers instability driven by a cylinder of non-circular cross section. 
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FIGURES 



FIG. 1. Time evolution of the asynchronous CML (|3|) with r = 3.0, e = 0.58 and system size 
L = 128. The turbulent sites with u < 1 are black while the laminar sites with u > 1 are white. 
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FIG. 2. Phase diagrams of one-dimensional CML for asynchronous updating (crosses) and 
synchronous updating (boxes). 
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FIG. 3. Order parameter m{t) versus time t for system sizes L between 2^ and 2^^ at r = 2.2. 
As expected the data falls on one line as long as t is much smaller than the escape time. 
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FIG. 4. Rescaling of the order parameter at r = 2.2 according to @ to obtain estimates for 
P/i'± and for z. Data for for system sizes L between 2^ and 2^^ collapses on one curve. 
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FIG. 5. Spatial correlation function Ci{t) for various times t at r = 2.2. Ci{t) approaches a 
straight line with slope 1 — r/ for large times, indicating an algebraic decay. 
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FIG. 6. Rescaled spatial correlation function at r = 2.2. The data for various times collapses 
on one curve if the exponents z and r] take the values in table ||. 
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TABLES 

TABLE I. Direct measurements of the critical exponents for system (^). The critical values Cc 
and the exponent z are found simultaneously by approaching the value of e where the scaling (|^ 
is best. The estimation of the other exponents are described in the text. The bottom line are the 
directed percolation exponents. 


T 

Cc 

z 

P 

A 

l/\\ 

ri 

2.2 

0.0195(2) 

1.58(2) 

0.28(1) 

0.16(1) 

1.51(2) 

2.6 

0.5096(2) 

1.59(2) 

0.28(1) 

0.16(1) 

1.49(2) 

3.0 

0.5870(3) 

1.60(3) 

0.28(1) 

0.15(2) 

1.50(2) 

DP 


1.57 

0.28 

0.16 

1.51 

TABLE II. 

The exponents and relations obtained from rescaling analysis using finite 

size seal- 

ing. The fourth column is an 

estimate of rj using the hyper-scaling relation (10). 


r 


z 

^ = 2i^ + l 

z 

V 

2.2 

0.26(1) 

1.57(1) 

1.52(2) 

1.56(2) 

1.53(2) 

2.6 

0.26(1) 

1.57(1) 

1.52(2) 

1.58(2) 

1.49(2) 

3.0 

0.25(2) 

1.58(1) 

1.50(3) 

1.58(2) 

1.53(2) 

DP 

0.26 

1.57 

1.51 

1.57 

1.51 
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